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Abstract 

We show an exact equivalence of the free energy of the g-state Potts antiferromagnet on 
a lattice A for the full temperature interval < T < oo and the free energy of the g-state 
Potts model on the dual lattice for a semi- infinite interval of complex temperatures (CT). 
This implies the existence of two quite different types of CT singularities: the generic kind, 
which does not obey universality or various scaling relations, and a special kind which does 
obey such properties and encodes information of direct physical relevance. We apply this 
observation to characterize CT properties of the Potts model on several lattices, to rule out 
two existing conjectures, and to determine the critical value of q above which the Potts 
antiferromagnet on the diced lattice has no phase transition. 
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The study of statistical mechanical models with magnetic field 0, temperature 
or both [|| generalized from real to complex values has yielded interesting insights into the 
properties of these models. Complex-temperature (CT) singularities have been shown to 
have, in general, rather different properties than physical critical points of spin models, 
including violation of exponent relations such as 7 7^ 7' || and violation of universality, as 
evidenced by lattice dependence of exponents (shown using exact results) 0. It is therefore 
of fundamental interest to understand better how CT properties of spin models are related 
to physical properties. Here we shall show an exact equivalence of the free energy of the 
g-state Potts antiferromagnet on a lattice A for the full temperature interval < T < 00 
and the free energy of the g-state Potts model (gPM) on the dual lattice for the interval 
—00 < exp(J/(k B T)) < — (g — 1), where J denotes the spin-spin coupling. This result 
leads to the important conclusion that there are two different types of CT singularities: 
the generic type, with violations of scaling and universality, and a special type, which is 
closely related to physical critical points and can be described by the same ideas of scaling, 
renormalization group (RG), universality classes determined as RG fixed points and related, 
for 2D models, to conformal field theories. Our observation also yields a way to determine 
some properties of certain theories with non-Gibbs measures, by relating them exactly to 
theories with Gibbs measures. This is a useful connection, since theories with non-Gibbs 
measures appear not just in the context of complex temperature or magnetic field, but 
also in a number of physical situations, such as quantum spin models in condensed matter, 
lattice quantum chromodynamics with finite chemical potential or topological charge, and 
the lattice formulation of the standard SU(2) x U(l) electroweak gauge theory. 

The g-state Potts model H-Pl with zero external field(s) on a lattice A is defined by 
the partition function Z = J2{a n } exp(— with the Hamiltonian 

n = -jJ2 * w (i) 

(nn 1 ) 

in standard notation, where a n = 1, ...,g, (3 = (/c^T) -1 , and {nn') denotes pairs of nearest- 
neighbor sites on A. We also define K = ft J, 

e K - 1 

a = z~ l = e K , x = (2) 

We consider d = 2 dimensional lattices, since in this case the model with spin-spin interac- 
tions along bonds maps to another with similar interactions on the bonds of the dual lattice. 
Physical applications of the 2D ferromagnetic (FM) Potts model for q = 3 and 4 include 
modelling properties of monolayers of gas molecules absorbed on substrates. The 2D Potts 
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ferromagnet is known to have a (7L q symmetry-breaking) phase transition which is second- 
order for 2 < q < 4 and first order for q > 5. The critical exponents and universality classes 
of the cases where the model has second-order transitions are well understood |J PH , but, 



aside from the q = 2 Ising case the free energy has never been calculated for general 
T. The Potts antiferromagnet has also been of interest because of its connection with graph 
colorings and the fact that, for certain lattices and values of q, it exhibits nonzero ground 
state entropy. 

Now consider the gPM on a 2D lattice A with iVo sites (0-cells), Ni bonds (1-cells) and 
N 2 faces (2-cells). The dual lattice = V(A) is defined by associating uniquely a (2— p)-cell 
of Arf with each p-cell of A. The partition function of the gPM satisfies the duality relation 

I BUI 

Z(A, q, a) = x^q N °- l -U^Z{k d , q, a d ) (3) 

where 

a + g — 1 1 

a d = V{a) = — , i.e. x d = - (4) 

a — 1 x 

It follows from ([3]) that in the thermodynamic limit, the singularities of the free energy of 

the gPM on A at a point a satisfy 

/(A ! Qi & ) sing — /(Ad j Qi ^d) sing (5) 

The key observation is that under duality, the complete physical temperature interval < 
T < oo, i.e., < a < 1 of the g-state Potts antiferromagnet (qPAF) on A is mapped to the 
CT interval — oo < a d < — (q — 1) on A^ and vice versa. (In contrast, the physical T interval 
of the g-state Potts ferromagnet on A, 1 < a < oo, maps under duality to another physical 
T interval, 1 < a d < oo, of the same model on A^.) 

From the above observation we can derive several interesting results. First, from eq. (|5|), 
it follows that there is an exact 1-1 correspondence between the CT singularities of the gPM 
on A d in the interval —oo<a d <—(q — l) and the physical singularities of the g-state Potts 
AF on A. In particular, if the gPAF on A has a phase transition at clpm-afm (here, PM = 
paramagnetic), then there is a corresponding CT singularity in f(A d ,q,a) at a = ai, where 

ai = V(apM-AFM) (6) 

and this is the leftmost CT singularity in the interval — oo < a d < (hence the label £). These 
transitions are of the same order. If the PM-AFM transition of the gPAF on A is continuous, 
then the singularity in / of the gPM at the CT point is the same when approached from 
the left or from the right and, writing f(A d ,q,a) sing ~ \a — ae\ 2 ~ ae , we have ag = a = a', 
where the exponent a = a' is the physical specific heat exponent of the g-state PAF at 
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the PM-AFM transition. (In general we cannot use duality to relate the f3 or 7 exponents 
for the staggered magnetization and susceptibility of the g-state PAF at olpm-afm on A 
and staggered or uniform (3 or 7 exponents at ai since the duality applies in the absence of 
external field(s).) Now for a model above its lower critical dimensionality, so that there is a 
symmetry-breaking phase transition, a general property of the CT phase diagram is that in 
the complex a plane, the continuous locus of points B where the free energy is nonanalytic 
(which determines the CT phase boundaries) is compact, i.e., the region sufficiently far from 
the origin is in the complex-temperature extension of the FM phase || [7|]. (Henceforth, we 
shall use (CTE)FM to denote the CT extension of the FM phase, and so forth for other 
phases.) It follows that if the gPAF on A has a PM-AFM transition, then since ag is the 
leftmost singularity in the a d plane of the g-state Potts model on A d , the CT phase to the 
left of ai is the (CTE)FM phase of the latter model. 

It may be helpful to illustrate these general results briefly for the known q = 2 case. For 
the square lattice, A = A^ and a sqA= 2,pM-AFM = v2 — 1, so that a S(jtq= 2/ = — (v^2 + 1). In 
the terminology of Ref. |7|], the CT phase to the right of a sqtq= 2/ is an O phase, meaning 
that it is not the CT extension of any physical phase. For the honeycomb (he) lattice, 
Q'hc,q=2,PM-AFM = 2 — v^, so that the leftmost CT singularity of the Ising model on the 
dual, triangular (t) lattice is at a t , q= 2,e = — v^3- For the Ising model on the honeycomb 
lattice, the phases on the left and right of ah c , q =2,e are the (CTE)FM and (CTE)AFM phases; 
for the Ising model on the triangular lattice, the phases to the left and right of at. q =2,e are 
the (CTE)FM and an O phase, respectively ||14|| . The analyticity of / on the he lattice in 
the interval — oo < a d < — 1 is now seen as being equivalent to the absence of a finite-T 
transition in the Ising AF on the triangular lattice. On the kagome lattice, ak ag , q =2,e = 
-3 1/4 (2 - V3)~ 1/2 [HI HI, separating the (CTE)FM and (CTE)PM phases on the left and 
right, and corresponding to the PM-AFM transition of the Ising AF on the dual diced lattice. 

A second type of behavior occurs if the gPAF on A has no finite-T phase transition but is 
critical at T = 0; this is associated with a part of the CT phase boundary B passing through 
a = in the CT phase diagram of this model. It then follows that for the g-state Potts model 
on A d a part of the respective B for that model passes through a d = —(q — 1). As examples, 
the Ising AF on the triangular lattice, and the q = 3 Potts AF on the square and kagome 
lattices all have zero-temperature critical points, so that there are respective singular points 
at the dual images, ah c , q =2,e = — 1, a SIM= 3^ = —2, and adi ce d,q=3,e = — 2 on the dual lattices. 
As discussed for the square lattice in Ref. ]17|], a S(M=3 ^ = —2 is the leftmost point where the 



CT phase boundary B crosses the negative real axis (see also Refs. |T(], |T] 

As the third type of behavior, if the g-state Potts AF on A has no finite-temperature 
phase transition and is also not critical at T = (which is manifested by B not passing 
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through a = in the complex a plane), then it follows that f(Ad, q, ad) is analytic in the CT 
interval — oo < < —(q — 1), which therefore must be part of the (CTE)FM phase of the 
qPM on Ad- This type of behavior occurs, for example, for the Potts model on the square 
lattice with q > 4. 

We consider now the q = 3 Potts model on the honeycomb and triangular lattices. 
Usually, in cases where there are no exact solutions, the locations of CT singular points 
are most accurately determined from series analyses; indeed, the complications due to these 
singularities were recognized in early work on series ||. However, for ag, one has a powerful 
alternative approach: to locate it using eq. @ and a precise determination of the physical 
PM-AFM transition of the gPAF on the dual lattice. This is our first application. For the 
q = 3 PAF on the triangular lattice, Monte Carlo and series analyses [0-0] led to the 
conclusion that this model has a weakly first order transition. The most recent study ||21|| , 
using Monte Carlo methods, has yielded a very precise determination of the transition point: 
T t = 0.62731 ± 0.00006, i.e., a PM -AFM,t = 0.20309 ± 0.00003. We infer that the q = 3 Potts 
model on the honeycomb lattice has a corresponding singularity at the CT point 

a hc , q=3/ = V(a PM _ AFM)t ) = -(2.76454 ± 0.00015) (7) 

i.e., Zhc, q =3,£ = a hc, q =3,e = ~ (0.36172 ± 0.00002) . This value is in agreement with, and is more 
accurate than, the value Zhc,q=3,£ = —(0.363 ±0.003) obtained from a recent low-temperature 
series analysis in Ref. |22| (where Zh c , q =3,e is denoted it_). If one accepts the reported 
result [Tj|-[^l|] that the PM-AFM transition of the q = 3 PAF on the triangular lattice is 
(weakly) first order, the same applies for the singularity at ah c , q =3,£, and, furthermore, these 
transitions have the same value of the latent heat. If one were formally to assign an exponent 
a to these transitions, it would thus be a = 1. Past experience shows that it can be difficult 
to distinguish a weakly first order transition from a continuous one, and indeed the series 
analysis in Ref. p2] found evidence for a continuous transition at z^ C)q= ^ with a ~ 0.5. 

A second application is to test a recent conjecture for an exact analytic value of the CT 
singularity at Zh c , q =3,i made in Ref. [p3fl , viz., that z hc , g=3 j = cos(2tt/9) - 3 1/2 sin(27r/9) = 



-.347296355.., i.e., ah c , q =3,e = —2.87938524... (arising as a root of the equation z 3 — 3z— 1 



0). Since the result of Ref. [^2 was 5a (where o = 0.003 was the uncertainty in the location 
of the singular point) away from this conjecture, it was concluded that this conjecture is 
unlikely but not impossible to be correct. We can strengthen this conclusion here: our 
determination of 27^,5=3^ decisively refutes the conjecture of Ref. [p3fl . 

As a third application, we use our recent study of the q = 3 Potts AF on the honeycomb 



lattice 24 1 , where we found that this model has no finite-temperature PM-AFM transition 



and is, indeed, disordered at T = 0, with nonzero ground state entropy So/ks = 0.5068(3). 



4 



From our discussion above, this implies that the free energy of the q = 3 Potts model on 
the triangular lattice is analytic in the interval — oo < a < —2. In passing, we note that 



this rules out yet another conjecture p5j , that there could be a singularity in this model 
at the largest negative root of the equation a 3 + 6a 2 + 3a — 1 = 0, viz., a = —5.411474... 
(associated with the completion of complex-conjugate branches of CT zeros labelled 6 in 
Ref. p5|). Our result also implies that for the q = 3 Potts model on the triangular lattice 
the leftmost component of the CT phase boundary B that crosses the real a axis must do 
so at a > —2. These results apply a fortiori to the case q > 4: in this range, the Potts 
AF on the honeycomb lattice has no finite- T phase transition and is disordered, with finite 
entropy, at T = 0, so that for each respective value of q, the corresponding Potts model on 
the triangular lattice is analytic in the range — oo < a < — (g — 1). These results are borne 
out by explicit calculations of CT zeros of the partition function for g-state Potts models on 
honeycomb and kagome lattices, to be presented elsewhere. 

A fourth application concerns the g-state PAF on the diced lattice (dual of kagome). Ref. 



[2| found CT singularities at z kag , q=3/ = -0.4023 ± 0.0005 and z kag , q=4/ = -0.42 ± 0.01 for 
the q = 3 and q = 4 Potts model on the kagome lattice. From the first of these, we calculate 
that the q = 3 PAF on the diced lattice has a PM-AFM phase transition at 

a dice d, q =3,PM-AFM = 0.1393 ± 0.0008 (8) 

For comparison, from the exact value for ak ag , q =2,i, we have adi C ed,g=2,PM~-AFM = 0.4354205... 
for the Ising case. Furthermore, we find the interesting result that the q = 4 Potts AF on 
the diced lattice has no finite-T phase transition and also is not critical at T = 0, since 

V{a kam =^) = -(0.18 ±0.02) (9) 

is negative. Remarkably, we have thus used information on a CT singular point of the Potts 
model on one lattice to derive a new physical result for this model on the dual lattice. From 
this we also determine the critical integer value of q for the g-state Potts antiferromagnet on 
the diced lattice above which it has no finite-temperature phase transition, viz., q = 3. 

In summary, we have used an exact duality property of a 2D g-state Potts model to relate 
certain complex-temperature properties of the model on a given lattice to physical properties 
of the Potts antiferromagnet on the dual lattice. This result gives useful information about 
the complex-temperature phase diagrams of 2D Potts models, and we have been able to 
use it to rule out two conjectures and to find the critical integer value of q for the Potts 
antiferromagnet on the diced lattice. 

This research was supported in part by the NSF grant PHY-93-09888. R.S. thanks Prof. 
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